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Abstract 

The flow over the zero-pressure-gradient So-Mellor 
convex curved wall is simulated using the Navier- 
Stokes equations. An inviscid effective outer wall 
shape, undocumented in the experiment, is obtained 
by using an adjoint optimization method with the 
desired pressure distribution on the inner wall as the 
cost function. Using this wall shape with a Navier- 
Stokes method, the abilities of various turbulence 
models to simulate the effects of curvature without 
the complicating factor of streamwise pressure gra- 
dient can be evaluated. The one-equation Spalart- 
Allmaras turbulence model overpredicts eddy viscos- 
ity, and its boundary layer profiles are too full. A 
curvature-corrected version of this model improves 
results, which are sensitive to the choice of a particu- 
lar constant. An explicit algebraic stress model does 
a reasonable job predicting this flow field. However, 
results can be slightly improved by modifying the as- 
sumption on anisotropy equilibrium in the model’s 
derivation. The resulting curvature-corrected ex- 
plicit algebraic stress model possesses no heuristic 
functions or additional constants. It lowers slightly 
the computed skin friction coefficient and the turbu- 
lent stress levels for this case (in better agreement 
with experiment), but the effect on computed veloc- 
ity profiles is very small. 

1 Introduction 

It has long been recognized that many turbulence 
models in use today are incapable of producing cor- 
rect physical behavior near curved surfaces. Much 
of this information comes from the use of boundary- 
layer codes (e.g., see Wilcox 1 ) applied to curved 
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flows with near-zero pressure gradient (e.g., So and 
Mellor, 2 Gillis and Johnston 3 ). The few applica- 
tions of Navier-Stokes codes to curved wall-bounded 
flows have generally been for cases with substan- 
tial pressure gradients, such as the U-duct test case 
of Monson and Seegmiller. 4 The use of a test case 
with pressure gradient can complicate the analysis 
by making it difficult to isolate the effects of curva- 
ture from the effects of pressure gradient. Also, in 
the Monson and Seegmiller case, one must contend 
with boundary layer separation and consequently a 
loss of two-dimensionality. 

The reason why Navier-Stokes codes have not 
been applied to zero-pressure-gradient cases (such as 
the So and Mellor case) is that in these experiments 
the outer wall shape has not been explicitly docu- 
mented. The only information recorded is that the 
outer wall was manually adjusted during the exper- 
iment to yield the desired (near-zero) pressure gra- 
dient along the inner curved wall. Additionally, in 
the case of So and Mellor, a local tangential jet was 
used at the outer wall near the start of the curve to 
maintain attached flow on this wall; while in the case 
of Gillis and Johnston, local boundary layer bleed- 
ing accomplished the same thing. From the stand- 
point of boundary layer methods, this information 
was adequate. However, these issues produce a sig- 
nificant challenge for the modeling of the entire two- 
wall setup, as is required in a Navier-Stokes simula- 
tion, and preclude a large amount of experimental 
data from being utilized. 

Today, with advances made in optimization meth- 
ods in CFD, it is now possible to find an outer wall 
shape that yields a specified pressure distribution 
on the inner wall. As a result, Navier-Stokes simu- 
lations can be relatively easily accomplished on test 
cases for which boundary layer codes were the only 
viable option in the past. 

In previous work, 5 three turbulence models were 
used to investigate the U-duct flow of Monson and 
Seegmiller. The three models employed were the 
one-equation Spalart-Allmaras (SA), 6 two-equation 
Menter shear-stress transport (SST), 7 and two- 
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equation explicit algebraic stress model (EASM). 2 * * * * * 8 
All models behaved similarly in the curved region, 
and all failed to predict the suppression of the tur- 
bulent shear stress caused by the convex curvature. 
Overall, the EASM was judged to be superior to the 
other two models for this flow field. However, as 
mentioned above, it is difficult to isolate the effects 
of curvature from other effects in the Monson and 
Seegmiller case. 

In the present study, we employ two of the above 
turbulence models (SA and EASM) to the So and 
Mellor case, which removes the complications of 
pressure gradient and boundary layer separation 
from consideration. Both of these turbulence mod- 
els include recently-developed curvature corrections, 
and can be run both with and without the correc- 
tions in place. We first describe an optimization 
method used to determine the outer wall shape, 
given the So and Mellor experimental inner wall 
pressure distribution. We then apply a Navier- 
Stokes code to the case. We attempt to answer the 
following questions regarding the isolated effect of 
curvature in zero-pressure-gradient flow: (1) how 

well do existing models without curvature correc- 
tion handle convexly-curved wall-bounded flow? (2) 
what aspects of the flow are missed, and how signif- 
icant are the missed effects? (3) how much improve- 
ment is gained by employing curvature correction 
terms to the turbulence models? 


3 Numerical Method and Turbulence Models 

The Navier-Stokes CFD code used in the current 
investigation is CFL3D, 12 a widely-used structured- 
grid upwind finite-volume method. Details about 
the code can be found in the User’s Manual refer- 
enced. 

The two turbulence models used are the one- 
equation Spalart-Allmaras (SA) model 6 and the two- 
equation explicit algebraic stress model (EASM). 5 
However, note that the EASM has an additional mi- 
nor modification, described in Gatski and Rumsey. 13 
Equation (4) in Ref. 5 is replaced by: 


7oT + tf 


(i) 


where 

7o = 7o — 1 (2) 

and 

7t =7i + l+( C ^ i _ < 7 1 )- (3) 

Gatski and Rumsey 13 showed that a source of er- 
ror in the EASM for curved flows was caused by the 
assumption of anisotropy equilibrium in the Carte- 
sian frame of reference in the derivation of the model 
directly from the full Reynolds stress model: 


Dijj 

Dt 


= 0 , 


(4) 


2 The Optimization Method 

The optimization method is described in Ander- 
son and Bonhaus 9 and Nielsen and Anderson. 10 
In summary, a fully discrete adjoint approach is 
used in an unstructured-grid framework to com- 

pute design sensitivities using either the Euler or 

the Navier-Stokes equations. The adjoint method 

includes the effects of the interior mesh sensitivities. 

A quasi-Newton optimization technique, referred to 

a s KSOPT, 11 is currently employed. 

In the adjoint approach for design optimization, 

a cost function is defined and augmented with the 

flow equations as constraints. In the present case, 
the cost function to be minimized is the pressure 
distribution on the inner wall, which is taken from 
experimental data. The shape of the outer wall is pa- 
rameterized with 28 design variables. To avoid hav- 

ing to contend with boundary layer separation along 

the outer wall, the optimization is conducted using 

the Euler equations and the method is run until the 

cost function reaches a suitable level of convergence. 


where = Tij/(2K) — S t j j 3 and I\ — r nn /2 is the 
turbulent kinetic energy. 

In the study of non-Newtonian constitutive rela- 
tions (e.g., Schunk and Scriven, 14 Souza Mendes et 
al. 15 ), a measure of relative rotation rate is based on 
the principal axes of the strain rate tensor. By as- 
suming a transformed form of Eq. (4) to hold in this 
principal axes frame, a new form of the EASM can 
be derived that takes into account the flow field cur- 
vature. This new form is termed EASM curvature- 
corrected (EASMCC). 

In the transformed coordinate frame, the following 
equation holds: 


where bij is the transformed anisotropy tensor. 
Written in the Cartesian frame, Eq. (5) becomes: 

= blk Q kj - Q lk b kj . (6) 

The tensor is related to the rate of rotation be- 
tween the principal axes (barred) system and the 
Cartesian (unbarred) system. 
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The method for implementation of EASMCC in 
2-D is as follows. The rotation rate tensor Wij = 
(dui/dxj — duj/dxi)/ 2 in the model is replaced 
by an “effective” Wjj = Wij - where the 

constant a 2 is defined by the pressure-strain cor- 
relation model. For the SSG model 16 used here, 
a 2 = ( 2 — C $)/2 and C 4 — 0.4. The tensor Hjj is 
given by 


Qij = 


0 Dot/Dt 
- DajDt 0 


(7) 


and DajDt is the Lagrangian derivative of the 
strain-tensor principal axes, given by 



tan 


-1 


y/5i 2 i + S? 2 - Sn 
S 12 


( 8 ) 


where Sij = (dui/dxj + duj /dx{)/2. This expression 
can be reduced to (see Ref. 17) 


Da 


Dt 2(Sl + S* 2 ) L 


c DSu c DSn 
an ~ tt ; ^12 ■ 


Dt 


Dt 


( 9 ) 


the model’s production term), which is not present 
in the EASMCC. 

In the current implementation of the SARC 
model, cm[1 — ft 2 ]Wv> a portion of the SA model’s 
production term, is replaced by Q,i[/ri — fn)Wu, 
where 

2 r * 

fr 1 = (1 + Cri )j— ~r [l - C r3 tan _ 1 (c r2 f)l - C rU 

(1 + r*) 

( 12 ) 

and c r i = 1, c r2 = 12. The constant c r 3 has been as- 
signed to be both 1.0 and 0.6 in Ref. 17, and Spalart 
and Shur admit that they are still experimenting 
with the function f r \. In the current study, we use 
both values, and show that c r 3 = 0.6 is the more 
appropriate choice for this case. The function r* is 
given by r* = S/W. The r term is computed using 


r = 2W ik S jk 


DSi 


Dt 


r + ( ffmn Sj n + tjmn.Sin'jQr 


/D\ 

(13) 

where the term represents the system rotation, 
and D — y / 0.5(5 2 + W 2 ). For 2-D flows and no 
system rotation, the expression for r reduces to 


In practice, a term is added to the denominator of 
Eq. (9) to avoid division by zero as well as to avoid 
spurious fluctuations in Da/Di in regions of very 
low gradient. Note that for a simple 2-D azimuthal 
flow with only a u$ component of velocity (a function 
of radius), the following relation can be derived: 

^ = \(S-W)sign(W l2 ), (10) 


where S = >/2 Sij Sij and W = y/2 WTjWij. This 
analytical function can be shown to hold in general 
only for H 2 very close to 1, where 1Z 2 is defined by 


U 2 


{W 2 } 

{S 2 } 


( 11 ) 


and {} represents the trace: {W 2 } = WijWji = 
-WijWij and {S 2 } = SijSji = Sij Sij. Eq. (10) 
has proved to be useful as a check (in regions where 
1Z 2 £3 1) on the more complicated numerics required 
to obtain DajDt exactly, but it is of limited use in 
general. We use the exact DajDt term given by 
Eq. (9) for all the results in this paper. 

A curvature correction for the SA model has been 
developed by Spalart and Shur, 1 ' and applied to a 
variety of flows in Shur et al. 18 This correction, 
Spalart-Allmaras for Rotation/Curvature (SARC) 
was similarly derived based on the rate of change 
of the principal axes of the strain rate tensor, but it 
also includes a heuristic function f r \ (that multiplies 


(Da\ 8VF 1S {S;, + sy 

m 1 (14) 

with DajDt given by Eq. (9). 

4 Results 

In the experiment of So and Mellor, 2 the curved- 
wall tunnel had an aspect ratio of 8 (depth of 48 in.) 
and the flow along the tunnel centerline was nomi- 
nally 2-D. Thus, 2-D computations are expected to 
adequately represent the flow field. The inner wall 
shape is defined by a series of 9 arc segments of vary- 
ing angle and radius. The initial radius of curvature 
is 10 in., and the final radius of curvature is 13.86 
in. The curved wall turns through a total of 150°. 
A detailed description of the inner wall shape can 
be found in So and Mellor, 2 The channel width 
is 6 in. at the inlet. An outer wall shape was ob- 
tained from the optimization program, which was 
run in Euler mode to obtain a shape such that the 
inner wall pressure distribution matched experiment 
throughout most of the curved region. A list of re- 
sulting outer wall points is given in Table 1. Inner 
wall points are given also, for reference. 

The grid employed in the Navier-Stokes computa- 
tions is shown in Fig. I. The grid size is 257 x 161, 
with a minimum normal spacing at the convex wall 
of 0.00015 in. This corresponds with a spacing in 
wall coordinates of approximately y + = 0.3. The 
grid extends from 24 in. upstream of the curved wall 
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to approximately 18 in. downstream of the end of 
curvature. The Reynolds number per inch is taken 
as 3.6417 x 10 4 , and the nominal Mach number at 
the inlet is M = 0.063. At the inflow boundary, the 
u- velocity profile is set based on the experimentally- 
measured skin friction and boundary layer thickness. 
The turbulence quantities are set to match the ex- 
perimental levels at the same location. At the out- 
flow boundary, pressure is specified at p/pref — 1 , 
and all other quantities are extrapolated from the 
interior of the grid. Additional details concerning 
the boundary condition specifications can be found 
in Rumsey et al. 5 

Slip-wall boundary conditions are applied at the 
outer wall in the CFD simulation. This boundary 
condition is consistent with the assumption used in 
the optimization method, and allows the simulation 
to be run without the complication of having to con- 
tend with tangential jet or bleed boundary condi- 
tions. At the inner wall, standard no-slip adiabatic 
solid wall boundary conditions are employed. 

For the remainder of the paper, we adopt a co- 
ordinate system with s measured along the inner 
wall in the flow direction and d measured normal 
to the inner wall. Thus, s represents the surface co- 
ordinates, or length measured along the inner wall 
surface. The boundary layer thickness at the inflow 
( s = 24 in.) is approximately 0.55 in., whereas at 
the start of the curvature (s = 48 in.) it is approxi- 
mately 0.95 in. Thus, at the start of curvature, the 
parameter S/R is roughly 0.095. According to Pa- 
tel and Sotiropoulos, 20 S/R < 0.01 represents very 
mild curvature, whereas 0.1 < S/R < 1 represents 
moderate to strong curvature. Therefore, the cur- 
vature for this case can probably be categorized as 
“moderate,” whereas the Monson and Seegmiller U- 
duct case (with S/R = 0.5) can be categorized as 
“strong.” 4 - 5 

In the experiment, the outer wall shape was set 
to allow a small pressure drop near the start of cur- 
vature, followed by a region of nearly constant pres- 
sure all the way to the end of the curved surface, at 
s = 79.43 in. Surface pressure coefficients are shown 
in Fig. 2 using two different turbulence models. The 
pressures match experiment very well over most of 
the inlet and curved wall segments. 

The effect of grid density on a typical solution is 
shown in Figs. 3 and 4. In these figures, the “fine” 
level is the full 257 x 161 grid, “medium” has ev- 
ery other grid point removed in both directions, and 
“coarse” has every other grid point removed again. 
For this flow, the skin friction shows about a 3 - 
5% difference between the coarse grid and fine grid 
solutions in the curved region, and roughly 1% dif- 


ference (or less) between results on the medium and 
fine grids. The turbulent shear stress shows a differ- 
ence between the coarse and medium levels, but no 
difference at all between the medium and fine lev- 
els. For all the remaining results in the paper, the 
medium level grid is employed. 

The effect of the SARC model constant c r 3 is 
shown in Fig. 5. Note that the cj levels are ref- 
erenced to the nominal velocity at the inlet, rather 
than the local “potential flow velocity at the wall,” 
as reported in Ref. 2. The experimental levels have 
been adjusted accordingly. The original SA model 
yields high cj levels over most of the curved- wall 
region. When SARC with c r 3 = 0.6 is used, cj 
levels agree well with experiment, but SARC with 
c r 3 =1.0 predicts cj levels that are too low. Unless 
otherwise noted, for all remaining SARC results, a 
value for the constant c r 3 = 0.6 is used. 

Surface skin friction results using all four versions 
of the turbulence models are shown in Fig. 6. EASM 
and EASMCC are both low near the beginning of 
curvature but are relatively close to experimental 
levels over much of the curved-wall region beyond 
s & 55 in.; EASMCC reduces the cj levels from that 
of EASM by only a modest amount . Overall, SARC, 
EASM, and EASMCC produce similar cj levels over 
most of the curved region in reasonable agreement 
with experiment. 

Velocity profiles in the bend are plotted in Fig. 7. 
In addition to results in the curved region, profiles 
are shown at the inlet ( s = 24 in. station), although 
experimental data is not available at this location. 
As mentioned earlier, the velocity profile is set at the 
inflow to match the experimental Cj , <5, and nominal 
velocity using law-of-the-wall relations. As seen in 
the figure, initial profiles at the inlet are essentially 
identical for all four models. In the curved region, 
the results begin to differ. The three models SARC, 
EASM, and EASMCC are very close to each other 
and are in good agreement with experiment, fknv- 
ever, the SA model predicts higher velocity levels 
over the first 20% of the boundary layer at all three 
stations. 

Turbulent shear stress profiles are plotted in Fig. 8 
for SA and SARC and in Fig. 9 for EASM and 
EASMCC. All shear and normal stress profiles, to 
be given below, are in the local body/normal coor- 
dinate system. Stresses in this frame are related to 
those in the Cartesian frame by the following rela- 
tions: 

u'u' = -(vVc — u'u' c )sin(2©) + uV c cos(20) (15) 

u'u' = u'i(' c cos 2 0 -f z/i/ c sin 2 0 + u'u' c sin(2©) (16) 
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t?V = v l y / c cos"0 -}- u* ti / c sin“0 — ■u / u / c sin(20) t (17) 

where the subscript c indicates Cartesian frame, and 
0 is the angle that the body tangent vector makes 
with the X-axis. In Fig. 8, the SA model significantly 
overpredicts the — uV levels in the curved region, 
whereas SARC agrees much better with experiment. 
The differences between EASM and EASMCC in 
Fig. 9 are much less marked. However, EASMCC 
is generally in better agreement with experiment, 
particularly for d/8 > 0.3, where the turbulence is 
suppressed to near-zero levels. 

The turbulent normal stresses are plotted for 
EASM and EASMCC in Figs. 10 and 11. Because 
EASM and EASMCC are nonlinear models, they can 
predict the normal stress differences between u f u f 
and v'v*. Results are in good agreement with exper- 
iment at the inflow and throughout the curvature 
region. The curvature correction in EASMCC has 
the effect of lowering the normal stress levels slightly 
from those of EASM. The u r u f and v*v* for SA and 
SARC are not shown. Linear eddy viscosity mod- 
els cannot predict the normal stress differences, al- 
though the ability to predict these differences is gen- 
erally not considered necessary for most thin shear 
flow applications. 

Finally, the velocity profiles are shown using wall 
coordinates in Figs. 12 - 14. The theoretical log-law 
curve plotted in these figures is due to Spalding. 21 In 
the experimental results of Fig. 12, it is noted that 
the effect of curvature is primarily in the “wake” re- 
gion beyond the log layer, where the u + levels are in- 
creased at successive stations downstream. The log 
layer itself remains unaffected by curvature. The 
SARC model overall reflects the correct trend, in- 
creasing in the wake region with downstream dis- 
tance in the curve. However, the effect is excessive 
when c r 3 = 1.0, and even the log layer itself is af- 
fected by the curvature correction and loses the cor- 
rect slope. When c r 3 = 0.6, a portion of the log layer 
retains the correct slope and only the region beyond 
« 100 is altered. The EASMCC model shows 
somewhat elevated wake levels of u + , similar in char- 
acter to the experiment, even with no curvature cor- 
rection. These levels are raised slightly through the 
use of the curvature correction in EASMCC. In both 
EASM and EASMCC, the log layer remains in good 
agreement with the theoretical slope. 

It is also instructive to return to check the original 
assumption made in the derivation of the EASMCC 
model. We already know from previous studies (e.g., 
Rumsey et al. 5 ) that the assumption Eq. (4) is not 
valid in regions of high curvature. We would now like 
to investigate the validity of the transformed equa- 
tion, Eq. (6). We do this by computing its actual 


value (the quantity Ujdb,k/dxj) at various locations 
in the converged solution, and comparing it to the 
quantity b ik Qkj - ttikbkj* 

Results are shown in Figs. 15 and 16 for Dba/Dt 
and Dbin/Dty respectively. The three successive 
curves in each figure represent results at the three 
stations in the curved region. It is shown in these 
figures that the assumption Eq. (6) is indeed valid 
in the curvature region, and is nearly exact in the 
lower part of the boundary layer. 

5 Conclusions 

This numerical study has yielded the following 
conclusions. The standard SA model (with no cur- 
vature correction) does not do a good job modeling 
the flow field with convex curvature. Eddy viscos- 
ity levels are significantly overpredicted, and velocity 
profiles are somewhat too full. The skin friction co- 
efficient in the curved region is overpredicted. The 
curvature correction in SARC significantly improves 
results, lowering eddy viscosity levels and bringing 
velocity profiles into better agreement with experi- 
ment. The best choice for the model constant c r 3 is 
0.6 for this test case. A value of c r 3 = 1.0 lowers the 
skin friction coefficient too much, and the log layers 
of the velocity profiles are significantly altered. 

The EASM and EASMCC models both do a 
reasonably good job predicting this flow field. 
EASMCC has an advantage over SARC in that there 
are no heuristic functions and no additional con- 
stants necessary. The model is derived by assum- 
ing anisotropy equilibrium in the reference frame 
defined by the principal axes of the strain rate ten- 
sor, rather than in the Cartesian frame for standard 
EASM. The modified assumption on the anisotropy 
tensor is shown to be valid in the curved region of 
the flow field. However, the resulting curvature cor- 
rection in EASMCC has only a minor effect for this 
case, slightly lowering the turbulent stress levels (in 
better agreement with experiment) and lowering the 
skin friction coefficient by a small amount. The ef- 
fect on computed velocity profiles is very small. 

Therefore it appears that some aspect of the 
EASM model enables it to perform reasonably well 
for this curved-flow case even without a curvature 
correction. Because the EASM is derived directly 
from the Reynolds stress model, it retains some 
of the invariance properties of the full differential 
form, even with the incorrect Dbij/Dt = 0 assump- 
tion. Thus EASM yields a better physical repre- 
sentation of the turbulence than lower-order models 
such as SA. By including the curvature correction 
(in EASMCC), all the frame-invariance properties 
are retained. This modification does improve cer- 
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tain details in this flow field, but overall the effects 
are relatively minor. 
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y, inches 



Figure 1. Grid for So-Mellor case, every other point 
of 257 x 161 grid shown. 
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20 40 60 80 
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Figure 2. Surface pressure coefficient, referenced to 
inlet conditions. 


Figure 3. Effect of grid size on surface skin friction 
coefficient, EASMCC model. 



Figure 4. Effect of grid size on turbulent shear stress 
profile at s = 71 in., EASMCC model (fine and 
medium curves are indistinguishable). 
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Figure 5. Effect of SARC model constant on sur- 
face skin friction coefficient, referenced to inlet con- 
ditions. 




Figure 7. Velocity profiles at $ = 24 in., s = 59 in., 
$ = 67 in., and $ = 71 in., referenced to inlet condi- 
tions (origin for each successive station is shifted 0.5 
units to the right). 



s, inches 


Figure 6. Surface skin friction coefficient, referenced 
to inlet conditions. 



0.000 0.002 0.004 0.006 0.008 
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Figure 8. Turbulent shear stress profiles for the SA 
and SARC models at s = 24 in., $ = 59 in., s = 67 
in., and s = 71 in., referenced to inlet conditions 
(origin for each successive station is shifted 0.002 
units to the right). 
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Figure 9. Turbulent shear stress profiles for the 
EASM and EASMCC models at s = 24 in., s = 59 
in., s = 67 in., and s = 71 in., referenced to in- 
let conditions (origin for each successive station is 
shifted 0.002 units to the right). 



0.000 0.005 0.010 0.015 0.020 
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Figure 10. u’u' turbulent normal stress profiles for 
the EASM and EASMCC models at s = 24 in., s — 
59 in., s = 67 in., and s = 71 in., referenced to 
inlet conditions (origin for each successive station is 
shifted 0.005 units to the right). 
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Figure 11. v ! v f turbulent normal stress profiles for 
the EASM and EASMCC models at s = 24 in., s = 
59 in., s = 67 in., and s = 71 in., referenced to 
inlet conditions (origin for each successive station is 
shifted 0.005 units to the right). 



Figure 12. Semi-log plot of the experimental velocity 
profiles. 
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Figure 13. Semi-log plot of the velocity profiles for 
SA and SARC (SARC with c r 3 = 0.6 is shifted 3 
units to the right, and SARC with c r 3 1.0 is 

shifted 6 units to the right). 



Figure 14. Semi-log plot of the velocity profiles for 
EASM and EASMCC (EASMCC is shifted 3 units 
to the right). 



D b n /Dt 

Figure 15. Comparison of actual Dbn/Dt with 
EASMCC assumption at s = 59 in., s = 67 in., 
and $ = 71 in. (origin for each successive station is 
shifted 0.004 units to the right). 



D b xz /Dt 

Figure 16. Comparison of actual Dbn/ Dt with 
EASMCC assumption at s = 59 in., s = 67 in., 
and s = 71 in. (origin for each successive station is 
shifted 0.004 units to the right). 
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Table 1. Wall points (in.) used for So-MeUor case 


xinner 

yinner 

xouter 

you ter 

.2400E+02 

.OOOOE+OO 

.2400E+02 

.6000E+01 

.3350E+ 02 

.OOOOE+OO 

.3406E+02 

.6000E+01 

.3580E4-02 

.OOOOE+OO 

.3668E+02 

.6004E+01 

.3 / 80E+02 

.0000E+00 

.3905E+02 

.5983E+01 

.3957E+02 

.OOOOE+OO 

.4119E+02 

.5969E+01 

.4116E+02 

.OOOOE+OO 

.4315E+02 

.5954E+01 

.4258E+02 

.OOOOE+OO 

.4495E+02 

.5958E+01 

.4388E+02 

.0000E+00 

.4661E+02 

.5978E+01 

.4506E+02 

.0000E+00 

.4815E+02 

.6171E+01 

.4615E+02 

.OOOOE+OO 

.4957E+02 

.6390E+01 

.4716E+ 02 

.1652E-02 

.5092E+02 

.6508E+0 1 

.4810E+02 

-.5451E-02 

.5220E+02 

.6466E+01 

.4898E+02 

-.6397E-01 

.5339E+02 

.6254E+01 

.4980E+02 

-.1668E+00 

.5448E+02 

.5904E+01 

.5057E+02 

-.3273E+00 

.5548E+02 

.5456E+01 

.5128E+02 

-.5392E+00 

.5639E+02 

.4941E+01 

-5194E+02 

-.7944E+00 

.5723E+02 

.4379E+01 

.5255E+02 

1081E+01 

.5799E+02 

.3784E+01 

.5313E+02 

1391E+0 1 

.5870E+02 

.3165E+01 

.5366E+02 

- . 172 1E+0 1 

.5935E+02 

.2530E+01 

.5416E+02 

-.2067E+01 

.5995E+02 

.1886E+01 

.5462E+02 

-.2428E+01 

.6051E+02 

.1236E+01 

.5506E+02 

-.2802E+01 

.6104E+02 

.5848E+00 

.5546E+02 

-.3189E+01 

.6155E+02 

-.6808E-01 

.5584E+02 

-.3588E+01 

.6202E+02 

-.7222E+00 

.5619E+02 

-.3999E+01 

.6248E+02 

-.1379E+01 

.5651E+02 

-.4419E+01 

.6290E+02 

-.2040E+01 

.5682E+02 

-.4849E+01 

.6331E+02 

-.2707E+01 

.5710E+02 

-.5288E+01 

.6369E+02 

-.3382E+01 

.5736E+02 

-.5737E+01 

.6404E+02 

-.4069E+01 

.5759E+02 

-.6195E+01 

.6435E+02 

-.4768E+01 

.5781E+02 

-.6663E+01 

.6464E+02 

-.5480E+01 

.5801E+02 

-.7142E+01 

.6489E+02 

-.6209E+01 

.5819E+02 

-.7631E+01 

.6510E+02 

-.6953E+01 

.5835E+02 

-.8133E+01 

.6528E4-02 

-.7713E+01 

.5848E+02 

-.8646E+01 

.6541E+02 

-.8490E+01 

.5860E+02 

-.9173E+01 

.6550E+02 

-.9282E+01 

-5869E+02 

-.9714E+01 

.6554E+02 

-.1009E+02 

.5876E+02 

-.1027E+02 

.6554E+02 

-.1091E+02 

.5880E+02 

-.1084E+02 

.6550E+02 

-.1175E+02 

.5881E+02 

-.1 143E+02 

.6542E+02 

-.1261E+02 

.5880E+02 

-.1204E+02 

.6530E+02 

-.1348E+02 

.5876E+02 

-.1266E+02 

.6514E+02 

-.1438E+02 

.5868E+02 

-.1331E+02 

.6495E+02 

-.1529E+02 

.5856E+02 

-.1398E+02 

.6471E+02 

-.1623E+02 

.5840E+02 

-.1467E+02 

.6443E+02 

- . 1 1 20E+02 

.5819E+02 

1538E+02 

.6411E+02 

-.1820E+02 

.5793E+02 

-.1612E+02 

.6373E+02 

-.1923E+02 

.5760E+02 

-.1G88E+02 

.6327E+02 

-.2028E+02 

.5719E+02 

-.1766E+02 

.6273E+02 

-.2136E+02 

.5671E+02 

1847E+02 

.6208E+02 

-.2246E+02 


xinner 

yinner 

xouter 

you ter 

.5613E+02 

-.1929E+02 

.6129E+02 

-.2356E+02 

.5544E+02 

-.2014E+02 

.6033E+02 

-.2464E+02 

.5462E+02 

-.2099E+02 

.5918E+02 

-.2567E+02 

.5363E+02 

-.2182E+02 

.5782E+02 

-.2663E+02 

.5246E+02 

-.2263E+02 

.5627E+02 

-2756E+02 

.5110E+02 

-.2345E+02 

.5454E+02 

-.2846E+02 

.4957E+02 

-.2433E+02 

.5263E+02 

• .2944E+02 

.4783E+02 

-.2533E+02 

.5055E+02 

-.3056E+02 

.4584E+02 

-.2647E+02 

.4826E+02 

-.3185E+02 

.4354E+02 

-.2780E+02 

.4572E+02 

-.3332E+02 

.4084E+02 

-.2935E+02 

.4288E+02 

-.3505E+02 

.3761E+02 

-.3122E+02 

.3962E+02 

-.3699E+02 
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